Correlations, spin dynamics, defects: the highly-frustrated Kagome bilayer 
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The SrCrgj,Gai2-9pOig and Ba2Sn2ZnGaio_7 P Cr7 P 022 compounds are two highly-frustrated mag- 
nets possessing a quasi-two-dimensional Kagome bilayer of spin | chromium ions with antiferromag- 
netic interactions. Their magnetic susceptibility was measured by local Nuclear Magnetic Resonance 
and non-local (SQUID) techniques, and their low-temperature spin dynamics by Muon Spin Res- 
onance. Consistent with the theoretical picture drawn for geometrically frustrated systems, the 
Kagome bilayer is shown here to exhibit: (i) short range spin-spin correlations down to a tempera- 
ture much lower than the Curie-Weiss temperature, no conventional long-range transition occurring; 
(ii) a Curie contribution to the susceptibility from paramagnetic defects generated by spin vacan- 
cies; (iii) low-temperature spin fluctuations, at least down to 30 mK, which are a trademark of a 
dynamical ground state. These properties point to a spin-liquid ground state, possibly built on 
Resonating Valence Bonds with unconfined spinons as the magnetic excitations. 

PACS numbers: 75.40.Gb, 75.50.Lk, 76.75.+i, 76.60.-k 
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I. INTRODUCTION 

A. Highly frustrated magnets 

Anderson's initial proposal of a Resonating Valence 
Bond (RVB) state was intended as a possible alterna- 
tive for the Neel ground state of a triangular network 
with Heisenberg spins coupled by an antifcrromagnetic 
interaction Q]. The RVB state, known also as a "spin- 
liquid" state because of the short-range magnetic corre- 
lations and spin fluctuations down to T = 0, was also 
proposed to explain the high-T c behavior of cuprates 
and, more recently, of the superconducting compound 
Na^CoC^ • 2/H2O |3|]. Although Anderson's conjecture 
was proven to be wrong for the triangular network 0, 
there is a growing consensus that the RVB state is the 
ground state of the so-called highly-frustrated networks 
[a, E| • Like the triangular network, the magnetic frustra- 
tion of these systems is exclusively driven by the trian- 
gular geometry of their lattice provided that spins are 
coupled through an antiferromagnetic (AFM) interac- 
tion. This is different from spin glasses, where frustration 
arises from the randomness of the magnetic interactions 
0- However, compared to a triangular network, these 
networks have a corner-sharing geometry, as in the two- 
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dimensional Kagome (corner sharing triangles, Fig. QJ, 
the three-dimensional pyrochlore (corner sharing tetra- 
hedras, Fig. QJ, and the quasi-two-dimensional Kagome 
bilayer (corner sharing triangles and tetrahedras, Fig. 0) 
lattices. The corner-sharing geometry introduces a "mag- 
netic flexibility" , thereby enhancing the frustration of 
these networks (thus the term "highly") and destabiliz- 
ing the Neel order. A highly-frustrated network, and 
the experimental counterpart known as highly-frustrated 
magnet (HFM), is therefore an ideal candidate to possess 
a RVB ground state. Moreover, since HFMs are insula- 
tors, the experimental investigation of their ground state, 
and the related theoretical modelling, is greatly simplified 
by the absence of phenomena such as superconductivity, 
charge ordering or itinerant magnetism. 

The original features of the ground state of the Heisen- 
berg two-dimensional (2D) Kagome network, which is the 
central issue of this paper, is already apparent through 
a classical description of their magnetism [3, Iii! LUj , 
yielding a dynamical ground state with infinite degener- 
acy. More precisely, the "order by disorder" mechanism 
selects a coplanar spin arrangement in the ground state. 
However, the system can swap from one coplanar spin 
configuration to the other by a rotation of a finite num- 
ber of spins, with no cost in energy. These zero-energy 
excitation modes, or soft modes, prevent a selection of 
one of the coplanar configurations [§3 and constitute a 
low-energy reservoir of magnetic excitations. It results 
in a finite entropy per spin in the ground state, without, 
still up to now, any definitive conclusions on a possible 
long-range order, even at T = 0. 

In a quantum description of the magnetism for Heisen- 
berg spins h on the Kagome network, the magnetic 
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FIG. 1: The Kagome (left) and the pyrochlore (right |T^|1 
corner-sharing lattices. The coordinance of each site is 4 and 
6, respectively. 




FIG. 2: Kagome bilayers in SCGO (left) and in BSZCGO 
(right). Two different sites are present in the Kagome bilayer 
lattice, one of coordinance 5, the other of coordinance 6. The 
magnetism of both compounds arises from the stacking of 
magnetically decoupled Kagome bilayers, ensured in SCGO 
by the presence of the non-magnetic singlets and in BZSCGO 
by a large inter-bilayer distance. 



designate by "disorder" , may release the frustration and 
stabilize an original ground state pel EHl I2I I2I I2I |25| . 
Experimentally, disorder is present in all HFMs and in- 
deed causes cooperative phase transitions at low tem- 
perature, reminiscent of spin-glasses or AFM systems. 
This is typically observed in the jarosites Kagome fam- 
ily, where the disorder is governed by spin anisotropy and 
by dipolar interactions |26( |95| . In the spin | Kagome 
[Cu3(titmb)2(OCOCH 3 ) 6 ]-H20 compound, the competi- 
tion between first- and second-neighbor interactions, re- 
spectively ferromagnetic and antiferromagnetic, yields an 
original ground state, with a double peak in the specific 
heat and plateaus in the magnetization under particular 
conditions |27|. In the S = 1 Kagome staircase M3V2O8, 
the competition between first and second neighbor inter- 
actions, combined to Dzyaloshinskii-Moriya interactions 
and spin anisotropy, produces a very rich field and tem- 
perature phase diagram 28j. In three dimensions (3D), 
dipolar couplings and interactions further than nearest 
neighbor yield long-range magnetic order at 1 K in the 
pyrochlore compound Gd2Ti 2 C>7 [2^. In Y2M02O7, lat- 
tice distortions have been proposed to relieve the frus- 
tration in the ground state [3(j]. However, some samples, 
like the pyrochlore Tb 2 _pYpTi 2 07 [3l|,|32| or the spinel 
ZnCr2C>4 |33|, display spin dynamics down to very low 
temperature. Anyway, the high spin values and/or the 
3D character of the frustrated geometry of these com- 
pounds remain far from the ideal Heisenberg Kagome 
case we are dealing with. Along with the spin-liquid 
ground state of the ideal systems, the plethoric num- 
ber of possible ground states in the presence of disorder 
(including the case of ferromagnetic interactions in the 
pyrochlore lattice, yiel ding the frustrated so-called "spin 
ice" ground state [34l l35j| , named so as it is equivalent 
to the order observed for hydrogen atoms in ice H 2 0), 
constitutes one of the rich aspects of the HFMs. 



ground state is disordered [Tlj. It is predicted to be 
RVB-like, built on a macroscopic number of singlet states 
[T3L llH . If any, a gap between the ground state and the 
first triplet state is expected to be fairly small, of the or- 
der of J/20 |k|j where J is the AFM interaction between 
nearest neighbor spins (94[. Most strikingly and similar 
to the classical case, the corner-sharing geometry and the 
half-integer spins lead to an exponential density of low- 
lying singlet excitations |16(, with energies smaller than 
the gap, and to a non-zero entropy per spin at T = 
[l7j . The same conclusions seem to emerge for the AFM 
pyrochlore lattice [l8j- The nature of the magnetic ex- 
citations, possibly unconfined spinons, is still an open 
question [l7 l fl9 | . 

As a major drawback, every additional contribu- 
tion to the nearest-neighbor AFM interaction, such as 
next-nearest-neighbor interactions, dipolar interactions, 
anisotropy, Dzyaloshinskii-Moriya interactions, lattice 
distortions or magnetic defects generated by the presence 
of magnetic impurities or spin vacancies, all of which we 



B. SCGO and BSZCGO 

Among all HFMs, the chromium-based spin | Kagome 
bilayer compounds SrCrg p Gai2-9pOi9 [SCGO(jp)1 |3r^| 
andtheBa 2 Sn2ZnGaio-7pCr 7 p0 2 2 [BSZCGO(p)] [33 are 
the most likely candidates to possess a spin-liquid ground 
state. This is not surprising in view of the low spin 
value, 2D character and low level of disorder of these 
compounds compared to other HFMs. They are in fact 
an experimental realization of a quasi-2D Kagome bilayer 
of spin I chromium ions with antiferromagnetic interac- 
tions (Fig. I2J). Moreover, the AFM interaction between 
nearest neighboring spins is dominated by a direct over- 
lap exchange resulting in a coupling of J ~ 40 K (Sec.lTlll. 
This coupling is almost two orders of magnitude larger 
than the typical disorder-related interactions, like single- 
ion anisotropy, estimated to 0.08 K [3^, |3^|, dipolar in- 
teractions ~ 0.1 K, or next nearest neighbor interactions 
< 1 K. These compounds, however, always have a small 
amount of substitutional disorder, since a Cr-coverage 




FIG. 3: Cr 3+ network of SCGO and BSZCGO with their 
oxygen (O) environments, along with the two Ga 3+ sites. 
The dotted lines represent some of the Ga-O-Cr hyperfine 
coupling paths. 



higher than p = 0.95 and p = 0.97 cannot be reached 
in SCGO(p) HI and in BSZCGO(p) 0, respectively. 
Hence there are always at least some percent of the mag- 
netic Cr 3+ ions which are substituted with non-magnetic 
Ga 3+ ions. These low amounts of spin vacancies, and 
the related magnetic defects that they produce, turn out 
to be insufficient to destroy the spin-liquid behavior of 
the ground state. We will show that local techniques like 
Nuclear Magnetic Resonance (NMR) and Muon Spin Re- 
laxation (/iSR) are the most suited to probe frustration- 
related properties in these systems, since they can bypass 
the magnetic contribution of these defects. 

Figure [3] presents a simplified chemical structure for 
the ideal p = 1 unit cell of SCGO and of BSZCGO, high- 
lighting the chromium Kagome bilayers |96j| . As it can be 
seen, BSZCGO is a pure Kagome bilayer, on the contrary 
of SCGO where there are two additional chromium sites, 
labelled Cr(c), between the bilayer. These pairs of Cr- 
spins are coupled by an AFM exchange constant of 216 K 
with a weak coupling to the Kagome bilayers (~ 1 K) 
|42|. They form an isolated singlet at low temperatures, 
and the resulting susceptibility is negligible at T < 50 K 
compared to the susceptibility of the Kagome bilayer in a 
pure compound. Moreover, the Ga/Cr substitutions are 
substoechiometric on these sites (see Sec. Ill C 3(1 . The 
full crystal structure is obtained by the stacking of the 
unit cells presented in Fig. |5J The magnetism of both 
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FIG. 4: Macroscopic static susceptibility Xmacro (T) of 
BSZCGO(0.97) measured with a SQUID under 100 G (left 
scale) and Xmacro (J) (right scale). Inset : ac susceptibility 
(78 Hz) below 2 K focusing on the spin-glass-like transition 
(T g ). 

compounds arises by magnetically decoupled Kagome bi- 
layers, ensured in SCGO by the presence of the non- 
magnetic singlets at low temperature (only weakly af- 
fected by the Ga/Cr substitution, see Sec. Ill C 3|) . and 
in BSZCGO by the large interbilayer distance of 9.4 A 
(in SCGO of 6.4 A). These compounds can therefore be 
considered as ideal two-dimensional systems. 

Figure 01 shows the typical macroscopic susceptibility 
Xmacro in these compounds. The linearity of Xmacro a ^ 
high temperature and the extrapolation of this line to 
XmLro = to a negative temperature is typical of AFM 
interactions and gives an order of magnitude of the Curie- 
Weiss temperature 9cw — zS(S + 1) J/3 (600 K and 
350 K for SCGO and BSZCGO). However, contrary to 
the case of "conventional" antiferromagnets ^3 1 no kink 
is evidenced in Xmacro around 0c w- This is actually one 
of the common signatures of the frustration in HFMs [j| : 
the magnetic correlation length cannot increase because 
of frustration, and mean field theory remains valid for 
T ~ 9c w- At low temperature, a spin-glass-like transi- 
tion occurs in both systems, around a freezing temper- 
ature T g ~ 3 K in SCGO(p) and 1.5 K in BSZCGO(p) 
(Fig. 0. 

However, unlike in other HFMs, the spin-glass transi- 
tion is unconventional. In fact, in conventional 3D spin 
glasses, the nonlinear susceptibility diverges at T g , the 
specific heat is proportional to T at low temperature, 
and T g is proportional to the number of defects. In the 
Kagome bilayers, however, (i) the nonlinear susceptibil- 
ity diverges at T g 01 > but T g increases when the number 
of spin vacancies decreases (Fig. QjJ , indicating that the 
freezing is related to an intrinsic property of the frus- 
trated bilayer. (ii) The specific heat is proportional to 
T 2 , as in AFM 2D long-range ordered systems. Its un- 
usual large value at low temperature unveils a high den- 
sity of low-lying excitations [321 13 1 an d its insensitivity 
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FIG. 5: Collection of spin-glass-like transition temperatures 
T g for SCGO and for BSZCGO versus Cr-concentration p (our 
data and from Ref. I^R l3Sl IBol IFlh . The line is a guide to the 
eye. 



to an external magnetic field suggests a large contribu- 
tion from singlet excitations, which is one of the most 
stri king features predicted for the S = -k Kagome lattice 
[45L l46t |47|. (iii) Magnetic fluctuations and short-range 
correlations, consistent with a 2D ma gne tic network, are 
encountered at low temperature [isl l49j. The Kagome 
bilayers therefore combine properties of 3D spin glasses, 
2D AFM ordering, 2D fluctuating magnetic states and 
original properties expected in S — \ Kagome systems. 

Given these two temperature scales of the macroscopic 
susceptibility, 9cw and T g , Ramirez proposed the defini- 
tion of a "frustration ratio", / = 0cw/T g , with / > 10 
in all HFMs H^. The Kagome bilayers SCGO and 
BSZCGO respectively display / - 150 and 230 [|l|, the 
largest ratios reported so far in compounds where the 
frustration is driven by corner sharing equilateral trian- 
gles Q . The large frustration, the Heisenberg spins along 
with the low disorder, make them the archetypes of HFM 
and maybe the best candidates for a RVB spin-liquid 
ground state. 

We present a review of the magnetic properties of 
SCGO(p) and BSZCGO(p) powder samples, which cov- 
ers a large range of Cr-concentration p. The compar- 
ison between samples of different concentration allows 
to identify the frustrated and disorder-related magnetic 
properties of the Kagome bilayer. Through local probes 
it is then possible to determine the true nature of the bi- 
layer 's static susceptibility (by NMR), as well as the bi- 
layer 's low-temperature dynamics and ground state mag- 
netic excitations (by /iSR). The outline of the paper is 
as follows. Section [H] is dedicated to the susceptibil- 
ity of the Kagome bilayer. It is shown that the sus- 
ceptibility of the bilayer of SCGO and of BSZCGO can 
be accessed directly through gallium NMR experiments 
E3 El HI HI Hi E3- Both susceptibilities exhibit a 
maximum in temperature at T ~ 40 K ~ J and decrease 
below this maximum down to at least J/4. This be- 



havior indicates that short-range magnetic correlations 
persist at least down to 10 K, consistent with the exis- 
tence of a small spin gap. The comparison of the Kagome 
bilayer 's susceptibility with x macro shows that the spin 
vacancies of the bilayer generate paramagnetic defects 
responsible for the low-temperature Curie upturn ob- 
served in the macroscopic susceptibility (both systems 
have also extra Curie contributions, coming from bro- 
ken Cr(c) spin-pairs in SCGO and from bond length dis- 
tribution in BSZCGO). Section ITTT1 is dedicated to the 
/iSR study of the low-temperature spin dynamics of these 

HFMs iii m mm. 

It is shown, that magnetic exci- 
tations persist down to at least 30 mK (the lowest tem- 
perature that could be accessed). This temperature sets 
an upper limit for the value of a spin gap. A qualitative 
and quantitative analysis of the data shows that these 
excitations are possibly coherent unconfined spinons of 
a RVB ground state. The energy scale T g would corre- 
spond then, in this phenomenological approach, to the 
signature of this coherent resonating state. A summary 
and concluding remarks can be found in Sec. I1VI 



II. FRUSTRATED VS. DISORDER-RELATED 
SUSCEPTIBILITY 

A. Gallium NMR in SCGO and BSZCGO 

The NMR experiments were performed on 69 Ga 
( 69 7/ 2tt = 10.219 MHz/T, 69 Q = 0.178 x 10~ 24 cm 2 ) 
and 71 Ga ( 71 7/2tt = 12.983 MHz/T, 71 Q = 0.112 x 
10~ 24 cm 2 ) nuclei in SCGO and BSZCGO using a §- 
t-7t spin echo sequence, where 7 and Q are, respectively, 
the gyromagnetic ratio and the quadrupolar moment of 
the nuclei. The gallium ions are located on two distinct 
crystallographic sites, labelled Ga(l) and Ga(2) in Fig.|3 
As pointed out in Ref. Hi 0, Hi Hi the interest of Ga 
NMR lies in the coupling between gallium nuclei with 
their neighboring magnetic Cr 3+ ions through a Ga-O- 
Cr hyperfine bridge (Fig-JSJ. In particular, the 69,71 Ga(l) 
nuclei are exclusively coupled to the Cr 3+ ions of the 
Kagome bilayer: to nine chromium ions of the upper and 
lower Kagome layers (labelled Cr(b) in Fig. and to 
three chromium ions in the linking site in between the 
two Kagome layers (labelled Cr(a) in Fig. [3J. Through 
the NMR of Ga(l), it is possible to probe locally the mag- 
netic properties of the Kagome bilayers of both HFMs. 
This is the central topic of the NMR study presented in 
this paper. 

Before going any further, we briefly recall some rela- 
tionship between the contribution of each gallium nucleus 
to the Ga(l) NMR line and its local magnetic environ- 
ment in order to underline what can be exactly probed 
through Ga(l) NMR. We suppose first that the suscep- 
tibility in the Kagome bilayer varies from chromium site 
to chromium site and label it, in a generic manner, by X- 
A Ga(l) at site i will contribute to the NMR spectrum 
at a position depending upon the number of the nearest 
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TABLE I: (Quadrupole parameters for the Ga sites in 
SCGO US E3 and in BSZCGO HJ. Notice that 69 vq = 
{ m Q/ 71 Q)v Q » 1.59 71 i/q. Quadrupole effects are therefore 
more pronounced for the 69 Ga isotope than for the 71 Ga. 





Ga site 


71 v Q (MHz) 


V 


SCGO 


Ga(l) 


2.9(2) 


0.005(6) 




Ga(2) 


20.5(3) 


0.050(35) 


BSZCGO 


Ga(l) 


3.5(1.0) 


0.60(15) a 




Ga(2) 


12.0(5) 


0.04(3) 



a Determined using point charge simulation. 



neighbors (NN) occupied chromium sites and their sus- 
ceptibility x- This corresponds to the shift K^> in the 
NMR spectrum for a gallium at site i 

occupied NN Cr(a,b) 

where A is the hyperfme constant (the chemical shift is 
neglected here). The average shift, K, of the NMR line, 
is simply related to the average K l over all the gallium 
sites. Hence, it is practically proportional to the average 
susceptibility XKag of the Kagome bilayer. On the other 
hand, the distribution of K l around K defines the mag- 
netic width of the Ga NMR line and reflects the existence 
of a distribution of \- As we show in the following sec- 
tions, the NMR width probes a susceptibility related to 
disorder in BSZCGO and in SCGO. 

Along with this hyperfme interaction, the 69 ' 71 Ga 
Hamiltonian in SCGO and in BSZCGO has a quadrupo- 
lar contribution. The quadrupole interaction of gallium 
nuclei in both SCGO and BSZCGO is a consequence of 
the coupling of the nucleus to the electric field gradient 
(EFG) produced by the surrounding electronic charges. 
Following the usual notations, the nuclear Hamiltonian 
may be expressed as 

H = - 7 Kl(I + if) H + ^[3J* -/*+,,(/*-/*)] , 

where Ho is the applied magnetic field, vq the quadrupo- 
lar frequency and < i] < 1 the asymmetry parameter. 
The principal axes of the shift tensor K are collinear with 
the direction of the nuclear spin operators I x , I y and I z . 

As shown in Fig.[3]the gallium sites of both HFMs have 
different crystallographic environments. Some simple ar- 
guments can be used to anticipate differences concerning 
their spectra (Tab.[IJ|: 

(i) The Ga(l) sites of both HFMs have a nearly- 
tetrahedral oxygen environment. Given this nearly-cubic 
symmetry, the EFG is small, i.e. vq is small. However, 
whereas the Ga(l) site of SCGO is only occupied with 
Ga 3+ ions, the Ga(l) site of BSZCGO is randomly occu- 
pied either by Zn 2+ or by Ga 3+ ions. As a consequence, 
point charge simulations of the EFG on the Ga(l) site of 
BSZCGO show that there is a large distribution of r\ and 
vq, contrary to the Ga(l) of SCGO. 




H(T) 



FIG. 6: Comparison of the 69 ' 71 Ga NMR spectra in SCGO 
and BSZCGO obtained by sweeping the field with a constant 
frequency m ~ 40 MHz, at T = 40 K. The arrows point the 
quadrupolar singularities for the Ga(l) and Ga(2), respec- 
tively located inside and outside the bilayers (Fig. ■ "CL" 
points the central line (i «-+ — i transition) singularities. The 
continuous and dashed lines show the shift of the in-range 
singularities of the Ga sites, from one system to the other. 
We cannot identify the first order satellites of the Ga(l) in 
BSZCGO on this spectrum. 



(ii) The Ga(2) site of SCGO is surrounded by a bi- 
pyramid of 5 oxygen ions, whereas the Ga(2) of BSZCGO 
lies in an oxygen-elongated tetrahedra. Both environ- 
ments have no cubic symmetry which yields a large 
quadrupolar frequency. For this site, point charge simu- 
lations yield a vq twice larger in SCGO than in BSZCGO. 

(iii) The ratio Ga(l):Ga(2) is 2:1 in SCGO, whereas 
it is 1:2 in BSZCGO. Consequently, the relative spectral 
weight of the Ga(l) and Ga(2) sites will be opposite in 
the two compounds. 

An extensive study of the 69 ' 71 Ga NMR spectra in 
SCGO was presented in Ref. Eol l53l It allowed to iden- 
tify the NMR lines and to extract the parameters of the 
nuclear Hamiltonian, among them the quadrupolar pa- 
rameters given in Tab. [I] It was shown that the Ga NMR 
spectrum of both isotopes is the sum of the Ga(l) and 
Ga(2) sites, plus a small extra contribution related to the 
presence of the non-stoichiometric gallium substituted on 
the chromium sites (labelled Ga(sub) in Fig.[SJ|. A similar 
analysis allowed to identify the NMR lines in BSZCGO, 
and to determine the related parameters of the nuclear 
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FIG. 7: Quadmpolar powder pattern simulations for both 
sites and isotopes with the parameters of Tab. [I] for frequen- 
cies vi = 40 MHz (left) and 84 MHz (right) in SCGO and 
BSZCGO. The area of each contribution is obtained consid- 
ering the stoichiometry of each site in the samples and the 
natural abundance of each isotope. The full NMR spectra 
are obtained by adding all the contributions (not shown). 



Hamiltonian |4l|, also listed in Tab.HJ 

Neglecting the Ga-substituted sites, the Ga NMR spec- 
trum displays four sets of lines corresponding to the two 
isotopes distributed on both Ga(l) and Ga(2) sites. In a 
powder sample, as used here, for a given site and a given 
isotope, the line shape results from the distribution of 
the angles between the magnetic field and the EFG prin- 
cipal axes. This yields singularities rather than well de- 
fined peaks (the so-called powder line shape), as in the 
Ga NMR spectra of SCGO(0.95) and of BSZCGO(0.97) 
presented in Fig. |SJ 

Powder-pattern simulations for two different r.f. fre- 
quencies (yi — 40 MHz and 84 MHz) are presented 
in Fig. [71 using the parameters of Tab. [I] They per- 
fectly agree with the experimental NMR spectra and 
clearly show that the BSZCGO spectrum is more intri- 
cate than the SCGO one, because of the more pronounced 
quadrupolar contribution. This is probably the reason 
why the Ga(sub) line cannot be resolved in BSZCGO, 
contrary to SCGO. 

Most of the SCGO data were acquired in a field-sweep 
spectrometer with a r.f. frequency v\ ss 40 MHz, be- 
cause then the 69 ' 71 Ga(l) NMR lines can be isolated, as 
shown in Fig. □ 40]. At high temperatures (Fig.Efr), the 
NMR line shape remains unaltered and shifts with de- 
creasing temperature. At low temperature (T < 50 K), 
the line broadens with decreasing temperature (Fig. 03). 
The 69 ' 71 Ga(l) shift is still qualitatively visible down 
to 5 K, where the Ga NMR signal is lost (see below). 
Fits of the low-temperature spectra (using a convolu- 
tion of a high-temperature unbroadened spectrum with 
a Gaussian) allow to extract quantitatively the shift and 
the width of the NMR lines for all the samples studied 
(0.72 < p < 0.95). All display similar shifts with temper- 



ature, but different broadening, the line shape increasing 
with Ga/Cr substitution at a given temperature [4fj|. 

In the case of BSZCGO, it can be seen from Fig. that 
it is advantageous to work at a frequency higher than 
40 MHz, in order to isolate the Ga(l) contribution |97| . 
Field-sweep 71 Ga NMR experiments were performed at 
a r.f. frequency vi « 84 MHz. At high temperature 
(Fig. [5^), a 71 Ga(l) shift is evidenced similar to the one 
in SCGO ;98]. At low temperature, the 71 Ga(l) and 
71 Ga(2) lines broaden and start to overlap. However, 
the two lines can be resolved by exploiting the different 
transverse relaxation times T 2 of the two gallium sites 
(T 2 , G a(i) ~ 20 /xs, T 2tGa{2) ~ 200 fis). The first step 
consists in using a large time separation between the two 
r.f. pulses (t = 200 [is) in order to eliminate the Ga(l) 
contribution. The isolated Ga(2) line is extracted using a 
gaussian convolution of the quadrupolar powder pattern. 
In the second step, a short time separation is employed 
between the pulses (r « 10 /is) so that the Ga(l) con- 
tribution is recovered. The Ga(l) line is finally isolated 
by subtracting the long r spectrum corrected in inten- 
sity to account for the shorter time separation employed 
[99| . Such a procedure can be employed down to 10 K. 
The width and the shift of the Ga(l) line can then be ex- 
tracted from a fit using a 71 vq = 3.5 MHz quadrupolar 
powder pattern convoluted by a gaussian line shape, and 
are little affected by the uncertainty on 71 vq. The ratio 
between the intensities of the Ga(2) and Ga(l) lines (in- 
tegrated area below the lines) was found to stay close to 
2 down to 10 K, in agreement with the 2:1 stoichiometric 
ratio for the Ga(2) and Ga(l) sites of BSZCGO. 

The temperature dependence of the intensity of the 
71 Ga(l) NMR line, which corresponds to the number 
of detected 71 Ga(l) nuclei, is presented in Fig. ^3 f° r 
SCGO. A similar trend is observed in BSZCGO, with 
larger error bars due to the preponderant Ga(2) contribu- 
tion. As temperature decreases the spin dynamics slows 
down, and results in a decreasing longitudinal relaxation 
time Ti (see Sec. IIII|I . hence in a progressive loss of the 
NMR signal below 15 K. Such a wipeout is reminiscent of 
the one observed in spin-glasses at temperatures near T g 
[5^ | . However, contrary to spin-glasses where ultimately 
the signal is recovered below T g , there is no evidence for 
such a recovery in SCGO, nor in BSZCGO. 



B. Susceptibility of the Kagome bilayer 

The Kagome bilayer susceptibility XKag is probed 
through the shift K of the 69 ' 71 Ga(l) NMR lines [H 
liH l54j . presented in Fig. ^] for the purest samples of 
both systems in a temperature range displaying no signif- 
icant loss of intensity (T > 10 K). The shift (and XKag) 
increases when the temperature decreases up to 45 K, 
where a maximum is reached, and then decreases again 
as the temperature is lowered. This behavior is common 
to both systems and is field- and dilution-independent 
[5^| . Most importantly, it is quite clear from Fig.llllthat 
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FIG. 8: 71 Ga spectra in SCGO (Hi = 2ttvi / 71 7) . a) High tem- 
perature. The 71 Ga(2) contribution appears as a flat back- 
ground in this field range, b) Low temperature, lines are 
broadened. The 71 Ga(2) contribution remains flat at the po- 
sition of the 71 Ga(l) line, which shifts to higher fields in this 
temperature region when the temperature decreases. 

there is a discrepancy between the low-temperature be- 
havior of XKag and of Xmacro- Instead of a maximum, the 
susceptibility probed by SQUID measurements follows a 
Curie-like law at low temperature. This e stab lishes that 
Xmacro is a two-component susceptibility 100], resulting 
from the sum of the Kagome bilayer's susceptibility and, 
as evidenced in Sec. Ill CI of a susceptibility related to 
magnetic defects (xdef) : 

Xmacro = XKag + Xdcf • (3) 

Concerning the high-temperature (T > 80 K) behav- 
ior of K, a phenomenological Curie-Weiss law, K = 
Cnmr/ {T + #nmr), is an accurate fit, yielding #nmr = 
440 ±5 K and 380 ± 10 K, respectively for SCGO(0.95) 
and BSZCGO(0.97). The "Curie- Weiss" constant Cnmr 
found from the fit is 20% larger in SCGO than in 
BSZCGO, which indicates that hyperfine couplings are 
stronger in SCGO, likely because of the shorter Ga-O-Cr 
bonds (Fig. |2J. Although very commonly employed to 
extract information from the HFMs susceptibility, this 
phenomenological law is a rather crude approximation, 
since the linear behavior of the inverse of the suscepti- 
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FIG. 9: 71 Ga spectra in BSZCGO. a) High temperature. The 
71 Ga(2) first order quadrupolar contribution appears as a flat 
background in this field range at the 71 Ga(l) line position, 
b) Low temperature, lines are broadened, o and □ are the 
short r and the rescaled long r spectra respectively. A is 
for the Ga(l) contribution, given by their subtraction. The 
dotted arrows point at the center of the Ga(l) line, which 
shifts to higher fields when T decreases. Continuous (dashed) 
lines are gaussian broadened 71 Ga(l) [ 71 Ga(2)] quadrupolar 
powder pattern simulation, with vq = 3.5 MHz (12 MHz) and 
■q = 0.6 (7? = 0.04). 



bility is expected to hold, within a mean field approach, 
only when T > 2#nmr, which is not the case here. From 
$nmr, one can therefore only grossly estimate the cou- 
pling constants J of the Kagome bilayers. In the ab- 
sence of any prediction for S = |, we used the high tem- 
perature series expansion of XKag, derived for a S = \ 
Kagome lattice, to correct #nmr by a factor 1.5 [(5(| ■ The 
coupling constant is finally extracted from the mean-field 
relation 6*nmr = l.5zS(S + 1) J/3. The coordinance is 
2 = 5.14 and corresponds to the average number of near- 
est neighbors for a chromium ion of the bilayer. The 
couplings then read J = 45 K and 40 K for SCGO and 
BSZCGO, respectively. These values are consistent with 
the cou plings o bserved in other chromium-based com- 
pounds [61], |62| , and they indicate that the Cr-Cr AFM 
interaction stems from the direct overlap between or- 
bitals of neighboring chromium ions, rather than from 
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FIG. 10: NMR signal intensity in several SCGO(p) samples. 
A loss of 69,71 Ga sites is observed for T < 5T 9 (~ 16 K). The 
line is a guide to the eye. 
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FIG. 11: K (left) and Xmacro (right) for SCGO(0.95). The 
NMR shift and the SQUID susceptibility do not follow exactly 
the same law at high temperature, since Xmacro also probes 
the susceptibility of the Cr(c)-Cr(c) spin pairs of SCGO (a 
detailed analysis can be found in Fig. 13 of Ref . l4ft) . 

an oxygen-mediated superexchange coupling. A direct 
interaction results in a coupling J very sensitive to the 
Cr-Cr distance (d), and follows a phenomenological law 

of 

SJ/Sd = 450 K/A . (4) 

The slightly stronger coupling found in SCGO compared 
to BSZCGO results then, following this viewpoint, from 
the shorter Cr-Cr bond lengths of its Kagome bilayer 
(Fig. . 

We now turn to the low-temperature behavior of XKag, 
and in particular we discuss the maximum around 45 K. 
The first possible interpretation of this maximum is the 
existence of a spin gap. Such a gap A is an important 
issue for the Kagome bilayer, since it would be the signa- 
ture of the existence of a singlet ground state, found in 
the theoretical quantum description of a S — \ Kagome 
layer ^{|. Its value is predicted to be A = J/20, and 



FIG. 12: Shift of the 71 Ga(l) line in SCGO and BSZCGO. 
Lines are GCC calculations (see text). Inset: Full scale down 
to 0. The line is the 36 spin S = \ cluster calculation men- 
tioned in the text (xed). 

should be of the same order for higher spins |101| . Since 
Ga NMR cannot access lower temperatures than 10 K, 
which corresponds for SCGO and BSZCGO to J/4, we 
cannot conclude whether the maximum in temperature is 
related to a gap. From an experimental point of view, its 
observation would require temperatures T < A in order 
to observe the exponential decrease predicted for the sus- 
ceptibility |47|. In the inset of Fig. ^1 we compare our 
experimental data to the susceptibility computed with 
exact diagonalization on a 36 spin S = | Kagome cluster 
[l9l | . We see that the position of the predicted maximum 
does not match the experimental one. Also, the sharp 
peak of the susceptibility of the calculation is not seen 
in XKag- Clearly, larger cluster sizes and the Kagome 
bilayer geometry are needed for a better comparison to 
the NMR data. For the time being, it may be safely 
concluded, that the behavior of the NMR shift, hence 
of XKag, is only consistent with a spin gap smaller than 
J/10. 

Rather than the signature of a gap, the maximum in 
XKag was assigned in Ref. |5J to the signature of a moder- 
ate increase of the spin-spin correlations of the Kagome 
bilayer which, however, must remain short ranged given 
the absence of any phase transition. This conclusion was 
confirmed by neutron measurements on SCGO for tem- 
peratures ranging from 200 K down to 1.5 K [63j, and by 
susceptibility calculations performed on the Kagome and 
the pyrochlore lattices with Heisenberg spins, in which 
the spin-spin correlation length is kept of the order of 
the lattice parameter These calculations, using the 
so-called "generalized constant coupling" method (GCC) 
|64|. consists in computing the susceptibility of isolated 
spin-clusters with zero-magnetic moments (triangles for 
the Kagome and tetrahedras for the pyrochlore) and in 
coupling these clusters following a mean field approach. 
In the case of the pyrochlore lattice, an exponential de- 
crease of the susceptibility is obtained at low tempera- 
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ture. In this case, the ground state of a cluster is non 
magnetic and the spin gap is between the S = ground 
state and the magnetic S = 1 excited states. In the case 
of the Kagome lattice, the ground state of the triangle is 
magnetic (Stotal = 2) wn i cn gives rise to the non-physical 
divergence in the susceptibility as T — > 0, due to the 
choice of a particular cluster j64| . The GCC simulations 
for the S = I Kagome and pyrochlore lattices, in between 
which the Kagome bilayer's susceptibility is expected to 
lie, are presented in Fig. ^| As shown, they agree with 
the data down to T = 40 K <C 0nmr, yielding a maxi- 
mum and a behavior around the maximum in agreement 
with the experimental findings. From the simulations, 
the values of C w are 260 K for SCGO(0.95) and 235 K 
for BSZCGO(0.97), which correspond to J = 40 K and 
J = 37 K, close to the previous values obtained through 
the high-temperature Curie- Weiss law. 

In conclusion, the GCC computation fits quite well the 
data down to T ~ J <C 6cw, where J w 40 K. It only 
shows that the strongest underlying assumption of this 
cluster mean-field approach is relevant, i.e., the spin-spin 
correlation length is of the order of the lattice parameter, 
which prevents any magnetic transition to a long range 
ordered state. 



C. Defect-related susceptibility 

In this section, we focus on the defects contribution to 
the susceptibility. The NMR linewidth and the SQUID 
data at low temperature show that the paramagnetic sus- 
ceptibility Xdef observed in Xmacro comes from the dilu- 
tion of the Kagome bilayer and from other kinds of de- 
fects, namely the Cr(c) pairs in SCGO and bond defects 
in BSZCGO. The Ga(l) NMR width not only enables to 
evidence that the vacancy of a spin on the network, i.e. 
the dilution, generates a paramagnetic defect, but also 
allows to shed light on the more fundamental question 
concerning the extended nature of the defects. We first 
present the experimental facts concerning the magnetic 
defects, as measured quantitatively through the NMR 
linewidth and the SQUID susceptibility Xmacro and then 
elaborate on their nature in both compounds. 
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FIG. 13: Ga(l) NMR width versus temperature for SCGO 
samples of different Cr-concentration. The width AH is nor- 
malized by the reference field Hi = vi/ 69 ' 71 'y to superimpose 
the results from the two gallium isotopes (as explicitly shown 
for p = 0.95). The solid lines are oc 1/T Curie-like fits. 




T(K) 

FIG. 14: Magnetic contribution to the NMR linewidth 
AH/Hi. For BSZCGO(0.97), the data for the 71 Ga(2) (□) 
are rescaled by a factor 6, corresponding to the ratio of the 
coupling constants, deduced from the high temperature shifts. 
• are for AH/H t [ 71 Ga{l)]. The lines are oc 1/T Curie-like fits. 
Inset: p-dependence of A# 71 [ 71 Ga(l)] at 50 K for SCGO(p) 
and BSZCGO (p). The lines are guides to the eye. 



1. NMR linewidth 

The Ga(l) NMR linewidth of SCGO and of BSZCGO 
was measured for various Cr-concentrations p. Figure ITSI 
presents the typical p-dependence of the linewidth in 
both samples — in the figure, only for SCGO. In agree- 
ment with the qualitative presentation of the raw spectra 
in Sec. Ill Al the width increases as temperature drops and 
is well described by 1/T Curie law (lines in Fig. I13|) . It 
is very sensitive to the lattice dilution, similarly to the 
low-temperature behavior of Xmacro (Fig. I16fl . but differs 
from the NMR shift. The perfect scaling of the widths of 
the two 69 < 71 Ga(l) isotopes (e.g. SCGO(0.95) in Fig.QS) 



underlines the magnetic origin of the low-temperature 
broadening. 

Figure [21 presents the linewidth obtained for 
SCGO(0.95) and for a comparable dilution in 
BSZCGO(0.97). Surprisingly, the Curie upturn in 
BSZCGO(0.97) is four times larger than in SCGO(0.95). 
This ratio cannot be explained by a higher hyperfine 
coupling constant for BSZCGO, since, as mentioned in 
Sec. Ill Bl the constants Cnmr obtained from the high- 
temperature Curie- Weiss analysis of the shift in SCGO 
and BSZCGO point to the opposite variation. Further 
insight concerning this mismatch between the two HFMs 
can be gained by plotting the widths as a function of 
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FIG. 15: 71 Ga(l) spectrum at 150 K in SCGO(0.95) {ui » 
40 MHz). The lines are quadrupolar powder pattern simu- 
lations with the parameters of Tab. [I] shifted with a value 
proportional to z', the number of NN Cr for a site. Their 
area is weighted by the probability P(z') when p — 0.95, pre- 
sented in the inset. 



dilution at a given temperature (inset of Fig. I14|) . In 
SCGO the width extrapolates to when p = 1, which 
shows that the width is only related to dilution. On the 
contrary, in BSZCGO the width reaches an asymptotic 
non-zero value for 0.86 < p < 1. As a first assumption, 
one could wonder whether the dilution of the lattice 
is larger than the nominal one. This can be ruled out 
since (i) the expected evolution of the line shape with 
p is found at 300 K (see the following) (ii) muon spin 
relaxation measurements indicate a regular evolution of 
the dynamical properties with p (Sec. IIII|) . Therefore, 
the increased low-temperature upturn in BSZCGO 
comes from dilution-independent paramagnetic defects, 
which are not present in SCGO. 

We may already rule out possible scenarios which could 
be responsible for the low-temperature broadening. For 
instance, the broadening generated by the suppression of 
chromium ions in the nuclear environment. Such a distri- 
bution is present in the Ga(l) NMR spectra of SCGO and 
of BSZCGO, but is only responsible for a minor broaden- 
ing of the line as shown in Fig.[JJfor SCGO(0.95). This 
broadening mechanism results from the fact that each 
Ga(l) nucleus has a probability P(z') (inset of Fig. IT5|) 
to have z' chromium neighbors (0 < z' < 12). As- 
suming that the Ga/Cr substitution does not affect the 
Ga-O-Cr hyperfrne couplings and that all the chromium 
ions have the same susceptibility, the shift of a Ga(l) 
nucleus surrounded by z' neighboring chromium ions is 
then proportional to z', according to Eq. ^ To con- 
struct the spectrum, each gallium nucleus is associated 
to a quadrupole line simulated with the parameters of 
Tab. [I] with a shift reflecting its Cr-environment and an 
intensity weighted by the total number of gallium nuclei 
having the same environment, i.e. weighted by P(z'). 
The simulated line perfectly matches the experimental 




20 40 60 80 
T(K) 

FIG. 16: Low-temperature macroscopic susceptibility 
(SQUID) for some SCGO Cr-concentrations. Inset: 
SCGO(0.95) data (open circles) fitted with Eq.El 

Ga(l) NMR line at 150 K (Fig. EH)). However, the broad- 
ening resulting from this distribution scales with the sus- 
ceptibility, hence with the shift K. Since K decreases 
at low temperature, this broadening mechanism cannot 
explain the low-temperature upturn of the width. For 
the same reason, a spatial distribution of the hyperfine 
constant, which yields also a width tx K, cannot justify 
the broadening observed and also has to be ruled out 
|l02j . Even if a such refined analysis is not possible in 
BSZCGO because of a stronger broadening than SCGO, 
similar conclusions are derived for this compound |56j . 
Hence, the explanation for the low-temperature broad- 
ening must be searched elsewhere (see Sec. Ill C 3l) . 



2. SQUID 

The SQUID measurements were carried out over a wide 
range of Cr-concentrations (0.29 < p < 0.97) for temper- 
atures down to 1.8 K (no difference was observed be- 
tween the Field Cooled and Zero Field Cooled suscep- 
tibility above the freezing temperature T g in a field of 
100 G). Figure [TBI presents the typical low-temperature 
macroscopic susceptibility of these HFMs. Like the NMR 
linewidth, the susceptibility exhibits a low-temperature 
upturn, increasing with growing dilution. In Sec. Ill Bl 
by a comparison between the NMR shift and the macro- 
scopic susceptibility, we established that Xmacro must be 
a two-component susceptibility, and, in particular, that 
it should possess, compared to the shift, a paramagnetic 
susceptibility Xdcf to explain its low-temperature upturn. 
The dilution-dependent upturn of x maC ro in Fig.llBlthcrc- 
fore establishes that Xdcf is related to the Ga/Cr substi- 
tution in SCGO(p). However, we also know from NMR 
that there also extra paramagnetic contributions, related 
to intrinsic defects in BSZCGO for example. 

In order to quantify the low-temperature dependence 
of the susceptibility, we follow Ref . |65| and fit x maC ro by 
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FIG. 17: Fits of the SQUID data (Eq.0. a) A come from 
Ref. IHol The lines are guide to the eye. b) The dashed line 
represents the Kagome bilayers defects contribution to Cdef (p) 
in SCGO(p), the other lines are guide to the eye. The grey 
area represent the Cr(c) contribution in Cd e f(p) for SCGO(p) 
(see text). Inset: same data as (a) for a broader range of p. 



a two component expression: 
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The first Curie- Weiss term roughly takes into account 
the contribution from the Kagome bilayer and governs 
the behavior of Xmacro at high temperature. The sec- 
ond term, the more relevant for this section, also a pri- 
ori Curie- Weiss, quantifies the contribution of Xdcf- A 
more accurate fit of Xmacro would consist in using a 
low-temperature susceptibility mimicking the NMR shift. 
Such a fit, however, yields the same qualitative results as 
presented below. 

From these fits (e.g. inset of Fig. Iltifl . an effective 
moment per Cr 3+ is extracted, p e ff = 4.1 ± 0.2/ in 
both SCGO ||3 and BSZCGO, close to the 3.87/^s ex- 
pected for S — | spins. In BSZCGO, the Curie- Weiss 
temperature is found to increase with p (Fig. I17fa ). in 
agreement with mean-field theory. From the linear vari- 
ation, following the pro cedure described in Sec. Ill Bl we 
extract J = 40 K [41| . consistent with our NMR re- 
sults. However, two linear regimes of 6(p) are observed 
for SCGO(p), below and above p ~ 0.55. This was first 
reported in Ref. |5(] and cannot be understood in terms 
of simple mean field theory. Taking into account a third 
term in Eq. corresponding to the susceptibility of the 
Cr(c) pairs in SCGO, derived in Ref. |4fJ and vanishingly 
small below 50 K, does not affect this result. 

The second term in Eq. is introduced to fit the low- 
temperature behavior of the defects. For all samples 
studied, it is close to a pure Curie law, within error bars. 
The values of the Curie constants Cdef extracted are pre- 
sented in Fig. 117b . A p-dependence qualitatively similar 
to the NMR width is found and Cd e f extrapolates again 
to for p = 1 in SCGO, whereas it displays a finite value 



in BSZCGO. When the Cr-concentration is p < 0.8, the 
dependency of Clef is qualitatively the same in both sys- 
tems (inset of Fig. 117b). To give an order of magnitude, 
this term represents, for the purest samples, a number of 
5 = 2 paramagnetic spins equal to the number of spin 
vacancies in SCGO(p) |4(| and to 15-20% of the number 
of Cr 3+ spins in BSZCGO(p) gj. 



3. Paramagnetic defects: discussion 

The NMR and the macroscopic susceptibility mea- 
surements indicate that there are paramagnetic defects 
in SCGO and BSZCGO. Some of them are related to 
the spin vacancies in the Kagome bilayers and affect the 
neighboring chromium ions in such a way that their over- 
all behavior is paramagnetic. These defects therefore 
yield information on how the correlated spin network of 
the Kagome bilayer responds to the presence of the va- 
cancy. In the following, we discuss this point as well as 
the presence of other types of defects in both compounds. 

The second type of defects result from lattice imperfec- 
tions in BSZCGO, and from the dilution of the spin-pair 
sites in SCGO. Unfortunately, in the case of BSZCGO 
the second type of defect dominate the paramagnetic re- 
sponse, and in this compound, on the contrary of SCGO, 
a fine analysis of the spin- vacancy defects cannot be car- 
ried out, even by Ga(l) NMR. 

a. Spin vacancies. In AFM systems such as the ID 
spin chains [MM! the quasi-2D spin ladders (6^ and the 
2D cuprates [69( it is now well established that a vacancy 
(or a magnetic impurity) generates a long range oscillat- 
ing magnetic perturbation and creates a paramagnetic- 
like, i.e. not necessarily cx 1/T, component in the macro- 
scopic susceptibility. The symmetric broadening of the 
NMR line observed in these systems is related to the 
oscillating character of the perturbation. The dilution 
effects of the Kagome bilayer would then be described by 
the same physics of these correlated systems. 

In SCGO, the symmetric feature of the Ga(l) NMR 
line indicates that spin vacancies induce a perturba- 
tion extended in space. This rules out descriptions of 
the defect in terms of uncorrelated paramagnetic cen- 
ters [6^, [?3 , since the broadening of the line would then 
be asymmetric. Hence, given the AFM interactions of 
the Kagome bilayer, the perturbation was assigned to a 
staggered polarization of the network [4(J ■ The exact di- 
agonalization of 36 spin i clusters on the Kagome lattice 
with spin vacancies by Dommange et al. confirms this in- 
terpretation |23| . They show that within an RVB ground 
state, the spin-spin correlations around a vacancy are en- 
hanced, and that the magnetization is staggered around 
it, a picture that should also hold for higher spin values. 
However, it seems that the image of a simple envelope 
of the staggered magnetization, like the Lorentzian now 
granted in High T c cuprates |69j | , is not correct here since 
this localization of singlets around the vacancies shifts 
the staggered cloud from the non-magnetic impurity. Al- 
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though the agreement is qualitative at the moment, fur- 
ther work is required within this framework for a quan- 
titative modelling of the NMR lineshape. 

b. Bond disorder in BSZCGO. We previously saw 
that in BSZCGO (p), low temperature macroscopic sus- 
ceptibility and NMR width data can be both satisfacto- 
rily accounted for, only if a novel p-independent defect- 
like contribution is taken into consideration, a result 
somehow surprising in view of the close similarity be- 
tween the Kagome bilayers of SCGO and BZSCCGO. 
The only major change lies in the 1:1 random occu- 
pancy of the Ga(l) site by Ga 3+ or Zn 2+ ions which 
induces different electrostatic interactions with the neigh- 
boring ions. As an example, distances to O 2- in a tetra- 
hedral environment vary from rg a 3+ -o 2 - — 0.47 A to 
rzn 2 +-o 2 - — 0-60 A which is at the origin of a change of 
the average Ga(l)-0 bonds, from 1.871 A for SCGO 
to 1.925 A in BSZCGO j^. Similarly, one expects that 
the Cr 3+ will be less repelled by Zn 2+ than by Ga 3+ , 
which certainly induces magnetic bond-disorder, i.e. a 
modulation of exchange interactions J between neigh- 
boring Cr 3+ . From Eq. ^ one can estimate that a modu- 
lation as low as 0.01 A in the Cr 3+ -Cr 3+ distances would 
yield a 10% modulation of J. 

The presence of non-perfect equilateral triangles clas- 
sically induces a paramagnetic component in the suscep- 
tibility of the frustrated units [7(J. The fact that we do 
not observe any extra Curie variation as compared to 
SCGO(p) in the average susceptibility, probed through 
the NMR shift K (Fig. EJ, indicates that such unbal- 
anced exchange interactions only induce a staggered re- 
sponse in the same manner as spin vacancies. Whether 
this could be connected with the localization of singlets in 
the vicinity of defects and a staggered cloud around them, 
like for spin vacancies, should be more deeply explored. 
This might highlight the relevance of a RVB approach to 
the physics of the Kagome network. Further insight into 
the exact nature of the defects likely requires a better 
determination of the bond disorder through structural 
studies at low T, to avoid the usual thermal fluctuations 
at room T . 

c. Broken spin-pairs in SCGO. The Ga/Cr substi- 
tution on a Cr(c) site of SCGO could in principle break 
a Cr(c)-Cr(c) spin pair and free a paramagnetic spin, 
which then could contribute to the paramagnetic upturn 
of Xmacro- However, this contribution turns out to be 
small. That the broken spin-pair susceptibility is negli- 
gible, can be deduced by comparing NMR and SQUID 
measurements for SCGO and BZSCGO. The starting 
point is to notice that the crossing between the varia- 
tion of the curie constants Cd e f of both samples deter- 
mined by the fits of Eq.03to the macroscopic susceptibil- 
ity fFig.H7b'). is not observed in the NMR width (inset of 
Fig. I14fl . Considering from our NMR analysis that there 
is a 20% difference between the hyperfine coupling, we 
can evaluate that the scaling factor between AH/H[(p) 
and C def (p) should be 20% larger in BSZCGO (p) than 
in SCGO(p). An error bar of 10% comes from the pos- 



sible difference between the bilayer susceptibility mea- 
sured with NMR and SQUID, since the Ga(l) nuclei 
do not pro be st oichiometrically the magnetic sites Cr(b) 
and Cr(a) |l03fl . We first determine the scaling factor 
[AH/Hi(p)]/C de{ (p) in BSZCGO (p), where 100% of the 
magnetic ions belong to the Kagome bilayers. With the 
same ratio, the AH/Hi(p) data for SCGO(p) yield there- 
fore a lower Cd e f(p) than the measured one, out of the 
error bars [dashed line in Fig. 117b], This difference is 
likely due to the Cr(c) pairs dilution in SCGO(p). Their 
contribution which, is not probed with the Ga(l) NMR 
[grey area in Fig.llTb] , corresponds to a O.OAp Cr 3+ S = § 
paramagnetic term, whereas a 2p/9 proportion would be 
expected in the case of a stoichiometrically substituted 
sample. This substoichiometric contribution of the Cr(c) 
site is consistent with neutron diffraction measurements 

0- 



III. SPIN DYNAMICS AND MAGNETIC 
EXCITATIONS 

A. Muon: an appropriate probe of dynamics 

We mentioned that an abrupt loss of the Ga(l)-NMR 
intensity occurs below 10 K (FigEIJ). Hence, Ga(l)- 
NMR cannot be employed to probe the low-temperature 
magnetic properties of SCGO and of BSZCGO. This 
is unfortunate, since the low-temperature spin dynam- 
ics reseal, in particular, clue information concerning the 
ground state of these systems. On the contrary, due to 
the smaller coupling of the muon to the electronic mo- 
ments and a shorter time window than the NMR one, the 
Muon Spin Relaxation (/xSR) technique has proven to be 
a front tool for the study of quantum dynamical states 
in a vast family of fluctuating systems. 

The electronic dynamics is commonly probed, either 
for NMR or for fj.SK, through the measurements of the 
longitudinal relaxation time, necessary to recover the 
thermodynamic equilibrium after the excitation of the 
nuclear spin or the muon spin systems. However, while 
the out-of-equilibrium state is reached using rf pulses 
in NMR, the muon spin system is already in an excited 
state. Indeed, muons (/i + , S = ^) are implanted in the 
sample, 100% spin polarized along the z axis. Therefore, 
the muon spins always depolarize to reach the Boltzmann 
distribution and fiSK can be performed in zero external 
field, contrary to conventional NMR [xj. The time de- 
pendence of their polarization P z (t) along the z axis is 
studied through their decay into a positron [72| and is di- 
rectly linked to both the spin fluctuations (the time corre- 
lation function of the local field is usually exponential, 
i.e., OMO)-H M (i))/(i^(0) 2 ) = exp(-z4)) and the (ran- 
dom) local field distribution, characterized by a width 
A/7 M (7^ is the muon gyromagnetic ratio). The longi- 
tudinal relaxation times is related to the spin-correlation 
function through 1/Ti - j* °°(S(0) • S(t)) cos(7 M # LF t)di. 
Zero field and longitudinal field (LF) //SR experiments, 
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where the external field Hlf is applied along the z axis, 
allow for instance to distinguish magnetic (randomly 
or ordered) frozen states from dynamical ones |73l l74j . 
Moreover, whereas specific heat is sensitive to all kinds 
of excitations, including low-lying singlets at low tem- 
perature in the Kagome bilayer samples |4fij |. the muons 
probe magnetic excitations only, with a specific range of 
frequencies (~ 10 9 Hz) sitting in a typical time window 
(10 ns-10 /zs) in between NMR and neutron experiments. 

We present here our /iSR stud y of the spin dynamics in 
SCGO(p) and BSZCGO(p) Hil58|. We first show that 
while conventional /iSR polarization functions can be 
used in the whole temperature range for the very diluted 
samples and in the high temperature regime (T > 3T g ) 
for the purest ones, they cannot be used satisfactorily in 
the low temperature regime in the latter case (p > 0.7). 
Moreover, in this first step, we qualitatively show that 
the magnetic state of all the samples is dynamical down 

We fur- 



to the experimental limit of 30 mK (Sec. IIII Bp 
ther use a model independent analysis of the data, simply 
taking the time necessary for the muon spin polariza- 
tion to decrease down to the value 1/e. The comparison 
between SCGO and BSZCGO clearly shows a correla- 
tion between the muon spin relaxation rate and T g , sug- 
gesting that this freezing temperature is not an impurity 
phase. On the contrary, it seems to be closely linked to 
the slowing down of the spin dynamics in the bulk sam- 
ple (Sec. IIII C|l . Finally, a phenomcnological model for 
the muon relaxation, based on sporadic dynamics due to 
spin excitations in a singlet sea, proposed by Uemura et 
al. is extended to all fields- and temperature- range. 
Its connection to the RVB picture is discussed, and we 
argue that such coherent states might mediate the inter- 
actions between impurities, which induce the spin glass 
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FIG. 18: Exponent (i{p) in BSZCGO(p) using the relaxation 
function P z (t) = exp[— (At)' 3 ] appropriate for the fast fluctu- 
ations limit at high temperature (T > 10 K.). 
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B. Conventional approaches 



High-temperature behavior 



FIG. 19: BSZCGO(0.97) data at 0.03 K. The lines are fits 
with different models (see text). We used A « 2 /zs , 4 fjs 
(y ~ 17 /is" 1 ) and 65 ^s _1 (y ~ 600 ^ts _1 )for the static, 
dynamic and "Uemura" cases respectively. 



At high temperature, a conventional paramagnetic be- 
havior is found for all the samples, with a stretched expo- 
nential variation of P z (t) = exp [-{\t) p }. FigurelTSlshows 
that as expected, (3 — > 1 in the dense magnetic cases jT^I 
(in our case with a high coverage of the Kagome bilayer 
lattice with Cr 3+ ) and /3 — ► 0.5 in the dilute cases [74J 
(i.e. with a low coverage). 

We can give here an estimate of the fluctuation fre- 
quency in this temperature range and of the NMR 
time window which would be required to measure this 
spin dynamics. At ~ 50 K, the muon relaxation rate 
is A - 0.03 (is- 1 and 0.01 /is" 1 in SCGO(0.95) and 
BSZCGO(0.97). In this appropriate fast fluctuating 
paramagnetic limit, the fluctuation rate can be estimated 
from v — ^/zJS/k B h — 2 x 10 13 s^ 1 49] using the cou- 
pling J w 40 K determined previously by NMR, and the 
average number z = 5.14 of Cr nearest neighbors. We ex- 



tract A = \J Xu/2 ~ 600 /us 1 and 300 fxs 1 , which cor- 
responds to an average field at the muon sites of 7000 G 
and 3500 G for SCGO and BSZSCGO respectively. 



Since the relationship between the longitudinal re- 
laxation time T\ measured in NMR and /zSR is 

T NMRJ T »SR _ [ A /(71 7 71 A )]2 j where 71 ^ „ 7Q k0 e//i S 

is the hyperfine coupling constant of the 71 Ga nuclei ob- 
tained from NMR experiments, we can estimate that 
rpNMR jrpjiSR _ 5 Thig ig consistent with our data 

^pNMR ^ 200 fjs at high temperature and establishes 
that below 5 K, the relaxation must be T^ MR ~ 1 /is, 
which is non-measurable as it lies out of the NMR time 
window. Such a small relaxation time is responsible for 
the wipeout of the Ga-NMR intensity (Fig. lTU|) . 
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2. Failure of conventional approaches at low temperature 
for the purest samples 

The derivation of zero field and longitudinal field /iSR 
relaxation functions is presented in Ref. |t3j IZ3 for text- 
book cases. We first rule out the possibility of some static 
freezing as the source of the muon relaxation. Two op- 
tions can be considered, (i) In randomly frozen (static) 
magnetic states, the major results are the following: first, 
the muon spin polarization displays a tail" in zero 
external field, i.e, P z (t > 5/ A) — > |. Indeed, -| of the 
frozen magnetic internal fields are statistically parallel to 
z and do not contribute to the muon spin depolarization. 
Second, P z (t) is "decoupled", i.e., does not relax, for lon- 
gitudinal fields Hlf > 10A/7 M . (ii) In the case of ordered 
magnetic states and powder samples, the " | tail" is still 
observed in zero field experiments for the same reasons, 
but oscillations appear in P z (t), which correspond to well 
defined local fields. 

One should notice that nuclear moments which appear 
static on the /iSR time window always contribute to the 
time-evolution of the polarization. A small longitudinal 
field of the order of a few 10 G (And ~ 0.1 /is -1 ) is com- 
monly applied to "decouple" this contribution. In other 
words, since the electronic and dipolar contributions are 
multiplied, such a contribution, if any, "disappears" and 
only the electronic contribution, not perturbed by such 
low fields, is measured. Hence, experiments performed in 
Hlf ~ 100 G are equivalent to a zero field measurement 
without nuclear dipoles. 

Figured shows P z (t) measured in BSZCGO(0.97) at 
0.03 K for H LF = 50 G and for H LF = 5000 G. Neither 
oscillations nor a ~ tail are observed when Hlf = 50 G. 
Although not emphasized in the figure, one can note that 
the initial polarization is Gaussian, which is even clearer 
in SCGO(p) jl^. I n a static case, such a shape is ob- 
served in dense magnetic randomly frozen systems, where 
the polarization is given by the Gaussian Kubo-Toyabe 
(KT) function We fitted our data at short times 

with this function in zero field, and plot the expected 
polarization when Hlf = 500 G in Fig. which is 
found to be completely decoupled. This contrasts with 
our experimental finding, i.e., the polarization at 5000 G 
is hardly decoupled. The lack of the ^ tail, the absence of 
oscillations and of decoupling effect show that the mag- 
netic state of BSZCGO(0.97) is fluctuating at 0.03 K. 
All these considerations remain valid in both SCGO(p) 
and BSZCGO(p) compounds when p > 0.6. 

The computation of the dynamical relaxation functions 
is based on a strong collision approximation 73] . Us- 
ing the corresponding Dynamical Kubo-Toyabe (DKT) 
function, analytically derived by Keren when v > A |75| , 
allows us to fit the zero field data on the whole time 
range (Fig. Again, the experimental data is much 

less sensitive to the applied field Hlf — 5000 G than the 
expected polarization. 

This so called "undecouplable Gaussian line shape" 
was first reported in Ref. for SCGO(0.89) and 



further in other Kagome compounds jH^, Eil |77| or 
spin singlet compounds like the doped Haldane chain 
Y2_ a; Ca a ;BaNi05 [T^. Uemura et al. proposed a relax- 
ation model |49j . presented further in this paper, which 
catches some of the facets of this relaxation. For now, we 
just notice that this model cannot justify our data for all 
fields (Fig- Ell ) even if it shows a lower decoupling effect, 
consistent with our (and previous) observations. Before 
going further with more complex models (Sec. 1111131) . we 
will first, in the following, give qualitative arguments to 
characterize the low temperature magnetism of our sam- 
ples. 

In order to illustrate the evolution of the properties 
with p, we present two typical low and high occupa- 
tions p of the Kagome bilayers in Fig. [201 which em- 
phasize the qualitative differences for various frustrated 
network coverage rates p. In addition to the high tem- 
perature behavior already commented, the evolution of 
P z {i) at low temperature is markedly different. For 
p = 0.97, the relaxation rate increases by more than two 
orders of magnitude to reach a temperature-independent 
value for T < T g m 1.5 K [Fig. 120Tb)]. with the un- 
decouplable character presented above. For a low cov- 
erage rate, p = 0.43, we also find a dynamical state 
but, at variance with the previous case, only a weak 
temperature-dependence is observed. Also, the polar- 
ization displays a square root exponential decay for all 
temperatures [Fig. l20T a1] and for any longitudinal field 
Hlf [Fig. l2Uf c1]. A weak plateau of the relaxation rate 
is still observed below 0.5 K but is no more present for 
p = 0.3, where P z (t) is found weakly temperature- and 
^-dependent . This is typical of the pure paramag- 
netic fast fluctuations limit for dilute magnetic systems, 
which contrasts with the results of SQUID measurements 
showing that strongly frustrated antiferromagnetic inter- 
actions are still present. This is the so called "cooperative 
paramagnetism" . 



C. Model independent basic analysis 

As a first step to a quantitative analysis, we estimate 
the muon spin relaxation rate A using the - point of the 
polarization P z (t) , as discussed in Ref. |57l This allows to 
single out information about two important issues using 
a simple model-independent analysis: the role of the spin 
glass-like transition and the impact of the spin vacancies 
on the dynamics. 

From high temperature, A increases by more than two 
orders of magnitude down to T g for both kagome bi- 
layers at their highest level of purity, to reach a relax- 
ation rate plateau At^o f° r T < T g 0, |^. Using a 
temperature-scale twice larger for SCGO(0.95) than in 
BSZCGO(0.97), the temperature dependence of A of both 
samples p erfec tly scales on the temperature axis as shown 
in Fig. [2] 104]. This ratio is very close to the 2.3(2) ratio 
between their freezing temperatures T g , which points to 
a link between the formation of the spin glass-like state 
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FIG. 20: a,b) temperature-dependence of the muon polariza- 
tion P g (t) in BSZCGO under weak H LF (10 G for p = 0.43 
and 100 G for p = 0.97). Notice the different time scales 
for both samples. c,d) //i^-dependence at base-temperature. 
The line for p — 0.43, ZF, is a square-root exponential times 
the Kubo-Toyabe function in zero field (And ~ 0.08 /^s -1 ), 
accounting for the nuclear dipole contribution. The other 
lines are fits described in the text. 



for SCGO(p < 0.89) in Ref.|57|. In a case of purely dipo- 
lar couplings of the muon with the electronic magnetic 
moments, the average of the local dipolar field is given by 
A cx (J^i 1 / r i ) oc p, where the sum is made over the Cr 3+ 
ions, rj is the distance between the muon and the ith posi- 
tion and (•) is the average over the muon positions. One 
therefore obtains X(p) ~ A.(p) oc p in the case of slow 
fluctuations and X(p) ~ 2A(p) 2 v / (v 2 + jhH^ F ) cx p 2 in 
the case of fast fluctuations, i.e., X(p) cx p v , with n < 2, 
which is not consistent with our data. This indicates 
that the relaxation is not induced by single spin excita- 
tions but rather collective excitation processes extending 
at least on triangles and/or hexagons. However these ex- 
citations still involve only a finite number of spins since 
the increase of spin vacancies only affect smoothly the 
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and the presence of fluctuations. This seems to be a quite 
common feature of various systems with a singlet ground 
state miH|7i|7l[7|. 

In Fig. [53 we report the variation of Xt^o(T) for vari- 
ous BSZCGO (p) samples and compare them to SCGO(p) 
[49T \Et\ . Although additional p- independent defects are 
dominant in BSZCGO (p), as evidenced by SQUID and 
NMR measurements, it is very surprising to find, for the 
two systems, a very similar quantitative low-temperature 
relaxation rate Xt^o with increasing Cr-concentration 
over the entire p-range studied. Therefore we conclude 
that only defects related to dilution of the frustrated 
magnetic network influence the relaxation rate. 

In a classical framework, a coplanar arrangement 
with zero-energy excitation modes, involving spins on 
hexagons, was proposed in the literature to be selected at 
low temperature [7!j . Inelastic neutron scattering exper- 
iments on SCGO and the spinel compound Zn&204 are 
indications in favor of such excitations |3.'l l48j . There- 
fore, one expects the muon spin relaxation to scale with 
the number of fully occupied hexagons, cx p 6 . In Fig. 1221 
the BSZCGO (p) data altogether with our SCGO(p > 
0.89) samples indicate that At^o(p) is well accounted for 
by adding a dominant p e term to the p 3 term proposed 



FIG. 21: Comparison of the evolution of the muon spin relax- 
ation rate A vs. T in BSZCGO(0.97) and SCGO(0.95), from 
a 1/e analysis. Inset: recovery of a part of the asymmetry 
at low temperature and long times in SCGO(0.95) below T g 
(log- log scale). 
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FIG. 23: Schematic representation of the creation and de- 
confinement of a spinon in a RVB ground state on a Kagome 
lattice. 

spin dynamics. On the contrary, a more pronounced ef- 
fect is observed in the S = | Kagome-like compound 
volborthite fll\ . 

D. RVB ground state with coherent unconfined 
spinons 

Finding a phenomenological model reproducing the po- 
larization P z (t) for all fields, all temperatures and all 
dilutions has been for long a pending challenge in the 
analysis of these fj.SH experiments in Kagome frustrated 
antifcrromagnets. In the quantum framework of a sin- 
glet ground state, one needs unpaired spins, i.e. isolated 
magnetic moments, in order to generate magnetic exci- 
tations responsible for the muon spin relaxation. Such 
excitations can be ascribed to unconfined spinons [l9| . 
for which the location of spin | varies in time with no 
loss of coherence of the excited state fFig. 1230. Hence, a 
given muon spin couples to a spin only a short fraction ft 
of the time t after implantation and one can use a model 
of "sporadic" field fluctuations to describe P z (t). 

This was the initial guess by Uemura et al. to account 
for the weakness of the field dependence of P z (t) at low 
temperature in SCGO(0.89) gjg, although the spins are 
S = §. In the absence of any theoretical predictions, we 
will just assume that the mechanism is comparable. This 
polarization is given by a "sporadic" dynamical Kubo- 
Toyabe function, P^ (ft, A, Hlf, v), which rewrites sim- 
ply Pf(t,fA,fH LF ,fv), A/7 P is the local field cre- 
ated on the muon site by a paramagnetic neighboring 
released spin with a fluctuation frequency v and a corre- 
sponding exponential time correlation function exp(— ut). 
A is therefore related to the muon location in the unit 
cell and its average is assumed to be p-independent for 
a given system but is expected to vary from BSZCGO 
to SCGO. The Gaussian at early times, the weakness of 
the field dependence and a dynamical relaxation down 



to are altogether related to the / factor and v ~ A. 
In addition to this sporadic relaxation, we found that a 
more conventional Markovian relaxation needs to be in- 
troduced to fit the long times tail (t > 3 /xs) for all fields 
and temperatures. We therefore write 

P z (t) = xPf(t, /A, f Hlf, fv) + (1 - x)e- yt , (6) 

where x represents the weight of the short time sporadic 
dynamical function, associated with the spinons dynam- 
ics. In the following, we first detail how all these param- 
eters can be reliably deduced from the data and sketch a 
physical picture consistent with our results. 

We first present our analysis for T <C T g in 
BSZCGO (p). In order to limit the number of free pa- 
rameters, we make the minimal assumption that the ex- 
ternal field does not influence the dynamics of the co- 
herent spinon term. The dynamics (v) and the average 
dipolar field created by a spinon on a muon site (A) are 
shared for all p and Hlf- For low fields, x is found close 
to 1 for the purest samples [Fig. I201 a.b)]. making A' a 
non-relevant fitting parameter. It also enables us to de- 
termine v and A when Hlf < 500 G, for various p. The 
parameter / is adjusted for each p and its variation ac- 
counts for the evolution of Ax— >o(p)- On the contrary, 
the high field data enable us to monitor the evolution of 
x with Hlf and to determine a value for A'. We could 
not extend the fits below p = 0.71 since the weak field 
dependence prevents an unambiguous determination of 
the parameters. 

We find perfect fits of our data [Fig. EJd)] with 
v ~ 1000 /is -1 , A ~ 350 /is -1 ~ 7 P x4 kG and an average 
value of / ~ 0.006. We find a nearly flat p-independent 
long time relaxation rate for T <C T g [A' ~ 0.05 ^s _1 , -j* 
in Fig. I24f d)]. An important finding is that x is of the 
order of p at low fields [Fig. l24f a.b)]. This may be related 
to the theoretical computations showing that the correla- 
tions are enhanced around spin vacancies in the Kagome 
lattice 23], which would destroy locally the spin liquid 
state. Besides, x decreases appreciably for Hlf ~ 10 kG 
whatever the value of p. We can associate this decrease 
to the existence of an energy scale ~ 1 K, which is of the 
order of T g . Finally, we observe a linear variation of / 
with p [Fig. l2"4T c)] and / tends to vanish around p ~ 0.5, 
a limit consistent with our classical approach fFig. I22fl . 
Indeed, since we find v ~ A when T — ► 0, the relation- 
ship f(p)A ~ \t^o(p) is expected. This indicates that 
even far from the substituted sites, the coherent state 
is somehow affected, e.g the density of spinons could be 
smaller. 

As suggested by the similar p- variation of Xt^o m both 
BSZCGO (p) and SCGO(p), we assume that the excita- 
tion modes are identical in both systems, i.e., / and v are 
kept the same for comparable p. We find A ~ 1200 /xs _1 
for our SCGO(p > 0.89) samples, in agreement with pre- 
vious work on SCGO(0.89) j43- It is quite rewarding to 
find that a common physical picture underlies all the sets 
of data at low temperature for both Kagome bilayers. 

We now extend the fits of the muon spin polarization 
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(open and full symbols). / is common for both systems. 



P z (t) to the whole temperature-range [Fig. I20f b)]. fix- 
ing / to its low-temperature value in order to limit the 
number of free parameters. It would also, phenomenolog- 
ically, mean that once a spinon is created, its deconfinc- 
ment process remains temperature independent, which is 
a quite reasonable assumption from a quantum point of 
view. As expected from the change of shape around T g 
[4^. , the weight x of the sporadic term decreases to 
finally enter a high-temperature regime [Fig. I24f b)] with 
an exponential muon relaxation (x — > 0). The similar- 
ity between the field- and temperature-dependence of x 
indicates that the sporadic regime is destroyed with an 
energy of the order of the freezing temperature T g . Fig- 
ure !22f b) displays a sharp crossover from one state to the 
other, between T g and 3T g , corresponding to the A steep 
decrease (Fig.^TJ, for both BSZCGO(p) and SCGO(p). 
For T > T g , v and A' [Fig. EUcQ] decrease by one order 
of magnitude up to 10T 9 . 

At high temperature, it is natural to think in terms of 
paramagnetic fluctuating spins. At lower temperatures, 
A' seems to diverge at T g and below, the weight of the 
exponential term at low fields (1 — x ~ 1 — p) could 
correspond to localized frozen spins, reflecting the glassy 
component measured by SQUID. It is noteworthy that T g 
does not increase with 1 — p but rather decreases, at the 
opposite of the case of canonical spin glasses. We could 
further confirm the existence of such a frozen component 
in SCGO(p = 0.95-0.89), since a clear recovery of a small 
part (~ 4%) of the asymmetry is found at long times [e.g. 
Fig. I21f inset)]. It is observed for t ~ 100 A -1 , which is 
out of the experimental range for BSZCGO (p). 

To summarize, the picture based on a coherent RVB 
state, which magnetic excitations are mobile fluctuating 
spins ^ on the Kagome lattice, explains well the data, 
provided that (i) these excitations can be generated even 



for T — > 0. This underlines the smallness of the "mag- 
netic" gap, if any, typically A < J/1000 ; (ii) an energy 
scale related to T g (fields of the order of 10 kG or, equiv- 
alently, temperatures of the order of T g , which vary very 
little with p) are high enough to destroy the coherent 
RVB type state ; (hi) the substitution defects are ac- 
counted for by an additional classical relaxation process. 



E. Spin-Glass-like transition 

The "intrinsic" spin-glass-like transition at T g is one of 
the most puzzling observation in these frustrated mag- 
nets. Indeed, the origin of such a spin glass state in a 
disorder-free system still awaits for a complete under- 
standing, although recent theoretical approaches catch 
some of the facets of this original ground state |80j . 
The most realistic model, presented by Ferrero et al. in 
Ref . 0, uses the so called "dimerized approach 0, 0] , 
i.e., considers a geometry of the Kagome lattice with 
two different kinds of equilateral triangles, correspond- 
ing to the Cr(b) layers in SCGO(p) and BSZCGO(p) 
(Fig. |3J). It predicts a spin glass- like transition in the 
S = \ Ka gome lattice, related to the freezing of the chi- 
rality, with T g ~ 0.5J'. Here, J' is the largest coupling 
in the Kagome planes (corresponding to the thick lines in 
Fig. |3J). Although 0.05 J' is of the order of magnitude of 
the experimental T g , the very similar values of the Cr-Cr 
bonds in both systems and Eq.0]do not allow to predict 
a factor 2 in T g . 

From our //SR results, we propose an other interpreta- 
tion which still awaits for theoretical confirmation. We 
find no energy gap for the spin | excitations. This could 
explain quite well why the transition to the spin glass 
state is fairly independent of p. Indeed, in this frame- 
work, spinons could mediate the interactions between 
magnetic defects localized around spin vacancies. This is 
corroborated by all the scaling properties in T/T g . The 
transition to the spin glass state would correspond then 
to the formation of the coherent singlet state rather than 
any interaction strength between defects. The bond de- 
fects inherent to BSZCGO (p) probably decrease the value 
of this coherent state energy scale since it is more favor- 
able to create localized singlets in this geometry. This 
is qualitatively consistent with the twice lower spin glass 
transition temperature in this family. 



F. Extension to other compounds? 

We already noticed that the /xSR undecouplable Gaus- 
sian line shape, altogether with a plateau of the relax- 
ation rate, have been reported in other compounds, and 
now elaborate about the possible links between these sys- 
tems. 

The other Ka gom e compounds, volborthite [Z]1IE2>E1 
and Cr jarosite [76j display the same properties as the 
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Kagome bilayers and our phenomenological model seems 
correct. 

Kojima et al. already noticed the similarity between 
SCGO and the doped Haldane chain Y2- x Ca K BaNi05 
and SCGO(0.89) ^1 & However, contrary to the 
Kagome bilayers, the relaxation function becomes more 
"conventional" when the system becomes purer, i.e., the 
muon spin relaxation function is a square root exponen- 
tial in the pure system at low temperature. It is not 
surprising since in the pure Haldane chain one expects a 
Valence Bond Crystal [T3. Therefore no magnetic exci- 
tation are allowed at low temperature and no depolariza- 
tion is expected in /iSR. However, mobile excitations are 
added when the chains are doped [85|. i.e. the muon spin 
relaxation becomes unconventional. This corresponds to 
the data and is consistent with our unconfined spinon 
phenomenological model. 

Finally, Fukaya et al. reported the same behavior in 
Sr(Cui_ x Zn l: ) 2 (B03)2 for x = and 0.02 86]. This sys- 
tem displays, theoretically and experimentally, an exact 
singlet dimer state and its intrinsic susceptibility van- 
ishes for T < 3 K HHH. Nonetheless, a CW term, 
corresponding to 0.72% of S = 5 impurities with re- 
spect to the Cu sites, is observed below 4 K [§3, which 
is finally comparable, from the macroscopic susceptibility 
aspect, to the former cases of pure volborthite and pure 
Y2CaBaNi05. However, according to us, there is today 
no clue about the possibility of mobile excitations in this 
magnetic network with a few percent of defects. 

We conclude that all these systems display (i) a theo- 
retical and/or experimental singlet ground state; (ii) a 
defect term observed through a low-temperature CW 
like upturn in the macroscopic susceptibility; (hi) a spin 
glass-like transition of this defect term at a temperature 
T g (but Sr(Cui_ a: Zn x )2(B03)2) GUI; (iv) a plateau of 
the muon spin relaxation rate below T g . 

It is important to notice that the clear experimen- 
tal correlation between the muon spin relaxation rate 
plateau and the spin glass-like transition is strongly 
against a pure muon-induced effect argument for this 
plateau. Even if the muon has indeed an effect on the 
physics of these compounds, the A plateau below T g 
shows that a bulk transition occurs at this temperature, 
i.e. that this freezing cannot be due to isolated impuri- 
ties. 

Theoretical approaches are now required to link a pos- 
sible unconhncd spinons state to a spin-glass-like behav- 
ior and a relaxation plateau of the implanted muon spins. 
They could therefore corroborate our phenomenological 
approach. 



IV. CONCLUSION 

With Hcisenberg spins and nearest neighbor couplings 
of the order of 40 K, the Kagome bilayers of SCGO and 
BSZCCGO are today the archetypes of Highly Frustrated 
Magnets. The comparison between SQUID and NMR ex- 



periments on the two slightly different compounds SCGO 
and BSZCGO allows us to isolate the intrinsic properties 
of the frustrated Kagome bilayer geometry and to un- 
derstand better the role of the (non)magnetic defects. 
We find that defects localized in the frustrated network 
induce a spatially extended response of the magnetic sys- 
tem, consistent with theoretical predictions. This re- 
sponse generates a Curie-like upturn in the macroscopic 
susceptibility at low temperature. However, local mea- 
surements of the susceptibility, using Ga NMR, show that 
the intrinsic susceptibility of the Kagome bilayers de- 
creases below 45 K, a temperature of the order of J. This 
behavior was also predicted theoretically and is consis- 
tent with short-ranged spin-spin correlations, on the or- 
der of the lattice parameter. We cannot give a definitive 
conclusion about the existence of a spin gap with NMR 
experiments because the gallium nuclei cannot be probed 
by NMR at low temperature. However a maximum value 
of ~ J 1 10 is found for such a gap. 

We measured the spin dynamics using Muon Spin Re- 
laxation experiments, down to 30 mK. We find that both 
systems SCGO and BSZCGO display a dynamical mag- 
netic state down to this temperature, unconventional for 
the purest samples. Very simple considerations show 
that the spin-glass-like state, measured at T g <C 8cw 
with macroscopic susceptibility techniques, is correlated 
to the slowing down of the Cr spin dynamics in the whole 
sample. The data are well accounted by a description 
based on unconfined spinons as the magnetic excitations 
and we suggest that this phenomenological approach re- 
mains valid for several spin singlet compounds at low 
temperature. Indeed, they all display striking similar 
aspects such as a spin glass transition of the magnetic 
defect channel in the macroscopic susceptibility and an 
undecouplable Gaussian line shape with a relaxation rate 
plateau of the muon polarization below this spin glass 
transition. 

In all these systems the fast depolarization of the muon 
down to the experimental limit of T ~ 30 mK underlines 
the existence of magnetic excitations at low temperature 
and hence the weakness of an hypothetical unconfined 
spinon gap. An other spin dynamics study with Neutron 
Spin Echo technique is in progress and could give more 
details about this T — > spin liquid state |89| ■ 
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quadrupolar frequency of the Ga(l) is smaller and the 
linewidth is proportional to the applied external field. 

[98] Notice that the high temperature analysis of both 
71 Ga(l) and 71 Ga(2) lines in BSZCGO allows us to com- 
pare their shifts, which yields an estimate of the chem- 
ical shift of 0.15(3)%. 

[99] We checked that the transverse relaxation time Ti is 
homogeneous over the whole 71 Ga(2) line. Its shape is 
hence used for the subtraction whereas its intensity is 
multiplied by ~ exp[2(r; — t s )/T2^]- 
[100] In the case of SCGO, Xmacro is a three-component 
susceptibility, the third contribution coming from the 
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[101] C. Lhuillier, Private communication. 

[102] A dipolar broadening due to the diluted paramagnetic 
defects |9l| would be to small to justify the broadening 
of the 69,71 Ga spectra. 

[103] Ga(l) nuclei are coupled to 9 Cr(b) and 3 Cr(a). The 
measured susceptibility is hence XKag.NMR oc 9xb + 3xa, 
where \b an d Xa are the Cr(b) and Cr(a) susceptibili- 
ties. On the other hand, the macroscopic susceptibility 
probes these susceptibilities with their stoichiometric ra- 
tio, i.e., XKag, macro oc 6xb + Xa E3 ■ A small difference 
between the ratios Xb/Xa would yield a slightly different 
ratio [AH/Hi(jp))/C del (p). 

[104] We attribute the factor 8 between A in both systems to 
a different coupling of the muon spin to the Cr spins. 
Transverse field experiments show for instance a larger 
linewidth in SCGO(0.95) than in BSZCGO(0.91) de- 
spite a lower defect term |4lll . On the other hand, A, 
which is mainly created by dipolar interactions in fiSH, 



is A oc Yl l/*" 3 ! where the sum is done over all the spins 
in the lattice and r is the distance between the muon site 
and these spins |73jl . Usually, the muon is located near 
an O 2- ion |92T| . In a very simplified approach, we com- 
puted A, considering that the muon would be located 
on the O 2- sites, as an average, and find A ~ 1100 /xs _1 
and 800 /xs" 1 in SCGO(l) and BSZCGO(l). This is in 
perfect agreement with the high temperature evalua- 
tion of A. Since A oc A 2 in the fast fluctuation limit and 
A oc A in the static limit, we expect, with this simple 
approach, A to be 1.5-2 times larger in SCGO(p). The 
missing factor 4 may either come from a more complex 
muon sites distribution, due to the O 2- ions around the 
Cr(c) sites, or to a different dynamics range as suggested 
by Neutron Spin Echo E^ . 
[105] According to us, no FC-ZFC magnetization data is pub- 
lished about Sr(Cui_ a: Zn l; )2(B03)2. 



